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Abstract 

We study the fractal geometry of 0(n) loop configurations in two dimensions by means of scaling 
and a Monte Carlo method, and compare the results with predictions based on the Coulomb gas 
technique. The Monte Carlo algorithm is applicable to models with noninteger n and uses local 
updates. Although these updates typically lead to nonlocal modifications of loop connectivities, 
the number of operations required per update is only of order one. The Monte Carlo algorithm is 
applied to the 0(n) model for several values of n, including noninteger ones. We thus determine 
scaling exponents that describe the fractal nature of 0(n) loops at criticality. The results of the 
numerical analysis agree with the theoretical predictions. 

PACS numbers: 64.60.Ak, 64.60.Fr, 64.60.Kw, 75.10.Hk 



I. INTRODUCTION 



The 0{n) model consists of n-component spins si = {s} , sf, . . . s^) on a lattice, with 
isotropic, i.e., 0{n) invariant couplings. The common form of the reduced Hamiltonian of 
the 0{n) spin model is 

H = -T^T.sl-s-^, (1) 

where the indices i and j represent lattice sites, and the sum is over all nearest neighbor 
pairs; J is the coupling constant, ks is Boltzmann's constant, and T is the temperature. 
Thus, the partition function of the model is 
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/ n exp s, ■ n ds~i (2) 
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where the spins are normalized such that This model includes as special cases the Ising, the 
XY and the Heisenberg model, for n = 1, 2 and 3, respectively. 

In the high-temperature limit, the bond weight w{si-Sj) reduces in first order to {1+xSi-Sj) 
with X = J/lksT). Thus the partition function takes the form 

^spin = / n (1 + ^■^'i ■ ■^1) n (3) 

The bond weight still satisfies the 0(n) symmetry implied by Eq. (^. According to the 
assumption of universality, the universality class of a phase transition is determined by 
only very few parameters including the symmetry of the spin-spin interactions. It is thus 
reasonable to expect (ij that the reduced Hamiltonian that corresponds with Eq. Q, namely 
H = — X]<jj>ln(l + xsi ■ Sj) with x = J/{kBT) not necessarily small, still belongs to the 
0(n) universality class in two dimensions. 

The 0(n) model (jHI) on the honeycomb lattice can be mapped on the 0{n) loop model 
P] on the same lattice, with a partition sum 

^loop = E ^''"n''^ (4) 

G 

where the graph G covers iV^ bonds of the lattice, and consists of Ni closed, non-intersecting 
loops. 

In the language of the 0{n) loop model, x is the weight of a bond visited by a loop, and 
n is the loop weight, and no longer restricted to be an integer. 



The research of 0(n) models is a subject of a considerable history, in which a prominent 
place is occupied by the exact results for the 0(n) loop model on the honeycomb lattice. 
These results include the critical points for —2 < n < 2, and the temperature and the 
magnetic exponent. 

Also the geometric description of fluctuations at and near criticality has a long history, 
which goes back to the formulation of phase transitions in terms of the droplet model Isf. 

n 

For the g-state Potts model (for a review, see Ref. y), it was found that the fractal dimension 
of Kasteleyn-Fortuin (KF) clusters |^ is equal to the magnetic scaling exponent yh- More 
generally, geometric Potts clusters can be defined by connecting neighboring, equal Potts 
spins by a bond percolationprocess. Several new critical exponents were found by Coulomb 
gas and other methods 0, Q, Isi . These exponents describe the geometric properties and 
the renormalization flow of this model. 

For the 0(n) loop model, the fractal dimension di of the critical loop gas has, to our 
knowledge, not yet been derived explicitly. However, a clue can be obtained from the 



relation lU| between the exponents describing random clusters of the tricritical Potts 

model, and those describing Potts clusters in the critical Potts model. From the equivalence 

fl 

of the critical 0{n) model and the tricritical q = n^-state Potts model one can thus 
associate the fractal dimension di of the critical 0(n) loop gas to the hull fractal dimension 

n 

of critical Potts clusters. The latter dimension was conjectured by Vanderzande |10l |. 

In the present paper, we focus on a fundamental non-thermodynamic scaling dimension 
behind the geometric properties of 0(n) loops which follows from the known correspondence 
with the Coulomb gas , and relate it to some exponents describing such properties. These 
exponents are exact and include the fractal dimension di of 0(n) loops at criticality. From 
a similar calculation, Duplantier and Saleur jl^ . llS^ derived the fractal dimension da of the 
interior of these loops. 

While theoretical predictions are available, thus far there is no numerical evidence in 
support of these, except for the Ising case n = 1 or g = 2 One of the reasons behind 
this situation may be that the 0{n) partition sum contains the number of loops. In the 
existing Monte Carlo algorithm for the loop model , the acceptance probability of a local 
update thus depends (for n 7^ 1) on the change of the number of loops due to the Monte 
Carlo move. The determination of this number is a nonlocal task and requires a number of 
operations that increases algebraically with the system size L. Critical slowing down can 



make this situation even worse, so that this way of simulation is restricted to rather small 
system sizes. 

Until now, a sufficiently efficient Monte Carlo algorithm for the 0(n) loop model has not 
been described. Therefore, in this work, we develop a new Monte Carlo algorithm, which is 
applicable to models with noninteger n > 1, and uses local updates. Although these updates 
typically lead to nonlocal modifications of the loop connectivities, the number of operations 
required per update is only of order one, and essentially independent of the system size. 

We then apply the algorithm to the critical 0{n) loop model and determine exponents 
of some geometric observables. The results agree with the theoretical predictions. 

The outline of the rest of this paper is as follows. In Sec. IHlwe show how a fundamental 
non-thermodynamic scaling dimension behind some geometric properties of 0{n) loops can 
be derived exactly from a mapping on the Coulomb gas, and how it relates to exponents 
describing some geometric observables. Section IIIII introduces the Monte Carlo algorithm. 
In Sec. II VI we apply the algorithm to the critical 0{n) loop model and determine exponents 
of some geometric observables. 



II. COULOMB GAS DERIVATION AND SCALING FORMULAS 



It is well known that geometric and fractal properties of 0{n) loops and various types 
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of critical clusters can be analyzed by means of a mapping on the Coulomb gas [l 
number of exact scaling dimensions was obtained by this technique, see e.g., Refs. 

Here we base ourselves on these analyses, which rely on a reformulation of correlation 
functions g{r) in the model of interest in terms of the Coulomb gas. The dimensions X{e, m) 
associated with such correlation functions are described by pairs of electric and magnetic 
charges, (cq, e^) and (mo,mr), separated by a distance r: 



X(e, m) 



(5) 



2g 2 

where g is the coupling constant of the Coulomb gas. For the critical 0{n) model it is given 

by 

1 n 

g = 1 -\ — arccos— (6) 



where we use a normalization that is in agreement with earlier literature but different from 



that used in Ref. 



16l (our g is four times smaller, and the charges differ by a factor two such 



that the X{e,m) are the same). 

Let us now consider the correlation function at criticahty describing the probabihty that 
two lattice edges separated by a distance r are part of the same loop. It decays with an 
exponent 2Xi where Xi is the 0(n) loop scaling dimension. The exponent Xi is described 
by a pair of magnetic charges mo = — = 1 and a pair of electric charges = = 1 — g. 
This leads to 

A-.^l-i (7) 

This dimension Xi is the renormalization exponent behind geometrical and fractal properties 
of 0(n) loops, just as the renormalization exponents Xt and X^. determine the thermody- 
namic singularities. 

In another application of the Coulomb gas technique we can explore corrections to scaling 
associated with the exponent X21 that describes the decay of the probability that two 
0{n) loops collide in two points separated by a distance r. The value of this exponent is 
determined by electric charges as above and magnetic charges itlq = —rrir = 2: 

X2i = l-- + ^ (8) 

This exponent becomes marginal at n = 2 and is thus expected to lead to poor convergence 
of finite-size data near n = 2. 

The physical relevance of Xi can be demonstrated by means of scaling arguments. The 
probability gi{r) that two points at a distance r lie on the same loop is, as given above, 
gi{r) ~ ar'"^^' . Let, at criticality, the fractal dimension of the loops be di. Thus, under 
a rescaling by a linear factor b, the length / of the loop decreases by a factor b'^', and its 
density increases by a factor b'^~'^i- . This determines the correlation in the rescaled system as 
giir /b) ~ a6^~^'^'r^^^' which is, as specified above, to be compared with gi{r/b) ~ afe^^'r"^^'. 
It thus follows that the fractal dimension of loops is 

di = 2-Xi (9) 

Let Pi{l) be the density of loops of length /. It is natural that, at criticality, Pi{l) depends 
algebraically on I, with an exponent denoted as pf. 

m)^p'. (10) 

Under a rescaling by a linear factor 6, the loop density is affected for two reasons: first, the 
loops decrease in length by a factor b'^'- = ; and second, the density increases by a factor 



6^ because the volume is reduced. Consistency requires that P/(/)(i/ = 6 ^Pi^lb"^ ^^)d{l}P' 
or Pi{l) = b-^'Piilb"^-^'). The requirement Pi{l) oc I'P' yields = (/6-2+x,)-k^-2+x,-2_ 
Matching the exponents shows that 

B = i + ^ (11) 

Since scaling implies that the divergence of the expectation value of the linear size of the 
largest loop g ) when the critical point is approached, and the 

actual length /max(a^) as expressed in lattice edges behaves as a power 2 — X; of the linear 
size, it follows that the largest loop length diverges as 

Ux(a;)oc(x,-a;)(^'-2)/^* (12) 

For L — s> oo and x > Xc, there exists an infinite spanning loop. Under a rescaling by a 
linear scale factor b its density increases by a factor b^^ while, as usual, the temperature 
field t oc X — Xc scales as t t' = The fraction si{x — Xc) of the edges covered by the 
spanning loop scales as si{W''{xc — x)) = b^^si{{xc — x). The choice b = {x^ — x)^^/^* leads 
to a constant on the left hand side of this equation, and after substitution in the right hand 
side, the scaling behavior follows as 

si{xc-x) (x{xc-x)^''y' (13) 

The finite-size dependence of the similar fraction sl of a system with finite size L at criticality 
can simply be found by rescaling the system to a given size, say 1. This leads to 

Sl oc L-^' (14) 

Including a correction to scaling, we may modify this into 

SL = aL-^'{l + bLy^ + ■■■) (15) 



where yi = 2 — X21 is a candidate for the leading correction exponent and a and b are 
unknown amplitudes. 

In analogy with magnetic systems, a susceptibility-like quantity xi can be defined on the 
basis of the distribution of the loop sizes as 

/max 

xi^T^mr (16) 

1=1 



According to the aforementioned scaling behavior, the largest loop in a critical system of 
finite size L has a length scale /max oc L"^"^' . Thus 

Xi{L) cx Y: (17) 
1=1 

Substitution of = 1 + 2/(2 — Xi) yields 

XiiL) ^ L'-'''' (18) 

Again, we can include corrections to scaling 

Xi{L)=cL'-''''{l + dLy' + ---) (19) 

with Ui as mentioned above, and c and d are unknown constants. 

Another correlation function of interest describes the probability that two sites on the 
dual triangular lattice separated by a distance r sit inside the same loop (i.e., not separated 
by any loop of the model). The exponent Xa describing the decay of this function at 
criticality was derived by Duplantier and Saleur as Xa = 1 — g/2 — 3/{8g). The fractal 
dimension of the interior of 0{n) loops is therefore da = 2 — Xa = l + g/2 + 3/{8g). The area 
inside a loop does not include the area inside loops enclosed by that loop. The exponent Xa 
thus also determines the finite-size scaling of the spanning loop. We are therefore interested 
in the fraction Sa of the dual lattice sites that sit inside the spanning loop. Scaling indicates 
that this quantity is subject to the following finite-size behavior: 

Sa{L) OC L-^'" (20) 

We furthermore define another susceptibility-like quantity Xa on the basis of the distribution 
Paid) of the area a (expressed in the number of enclosed sites of the dual lattice) of the 
interior of the loops as 

Xa{L)^L-'Y.Paia)a' (21) 

a 

which is expected to scale as 

Xa{L) OC L^-^^'' (22) 



III. MONTE CARLO ALGORITHM 

In the existing Monte Carlo algorithm for the loop model [14], the local updates require 
time-consuming nonlocal operations as explained above. To get rid of these we adopt the 
following procedure. 



FIG. 1: Representation of a loop configuration with Ising spins 



As a first step of such an algorithm for the simulation of the 0{n) loop model on the 
honeycomb lattice, we represent the loop configuration by means of Ising spins on the dual 
lattice, which is triangular. The loops are just the interfaces between neighboring spins of a 
different sign. We restrict ourselves to systems with periodic boundary conditions, so that 
the interfaces indeed form a system of closed, nonintersecting loops loops on the honeycomb 
lattice. This is illustrated in Fig. ^ using a loop configuration which consists of 2 loops 
and 16 bonds, shown as bold lines. This graph contributes a weight x^^n^ to the partition 
function. The loops can of course be represented by two opposite Ising spin configuration, 
but this degeneracy has no further consequences for our line of argument. 

We now show how one can update the loop configuration by means of local Metropolis- 
type updates of the Ising spins representing the loop configuration. It works only for n > 1. 
One step of importance sampling is realized by the following operations, which are to be 
repeated cyclically: 

1. For each loop, assign its color to be either 'green' with probability 1/n or 'red' with 
the remaining probability 1 — 1/n. 

2. Randomly select an Ising spin on the dual lattice. 

3. Check if the spin is adjacent to a red loop segment. If so, do nothing; if not, update 
the spin using the Metropolis probabilities, with Ising couplings K determined as 
e = X. 

4. Repeat steps 2 and 3 until the number of update attempts is equal to the number of 
sites of the dual lattice. 



5. Perform a sweep through the Ising system to find all loops. 

6. Repeat steps 1 to 5 a fixed number — 1 times. 



7. Sample the data of interest from the loop configuration. 

A Monte Carlo run consists of many of these cycles, each of which thus includes Us Metropolis 
sweeps, new random assignments of loop colors, and the data sampling procedure. Each of 
these steps satisfies the conditions of detailed balance, so that the algorithm should indeed 
generate configurations in accordance with Boltzmann statistics. Tests confirm that the 



simulation results agree with those of the existing algorithm 1J|. Since this algorithm 
assigns colors to the loops, we refer to it as "coloring algorithm". 



IV. SIMULATION 



With the help of the algorithm described above, we simulated the Oin) model on the 

n 

honeycomb lattice at the exactly known jjj critical points = (2 + (2 — n)^/^) in order 
to check the theoretical predictions described in Sec. |Tll 

We first applied our algorithm to the case n = 1.5, using 12 system sizes in the range 
8 < L < 128, where L is the linear size of the dual triangular lattice, and periodic boundary 
conditions. 

For each system size, a run was executed with a length of 4 x 10'^ Monte Carlo cycles 
after equilibration of the system. Each cycle included Ug = 5 sweeps and loop formation 
steps, and one sampling as described above. Statistical errors were estimated by means of 
binning in 2000 partial results. 

The data sampling included the density Pi{l) of loops of length /. The results for the 
system of size L = 128 are shown in Fig. |21 It displays a substantial interval of algebraic 
decay, as described by Eq. (fTUj) . 

The fractions s/ and Sa, and the susceptibility-like quantities xi and Xa were also sampled. 
The finite-size dependence of these quantities at criticality is shown in Table HI 

These quantities are well described by power laws as a function of the lattice size L for 
sufficiently large L. This is just as expected on the basis of finite-size scaling as expressed 
by Eqs. (HI), (UHl), (EOl) and ((221). This behavior is illustrated in Figs. El H El and El 

n 



TABLE I: numerical data for si{L),xiiL) and Sa{L),Xa{L) for different system sizes L at the 
critical point of 0(n = 1.5) model 



L 


Sl 


XI 


Sa 


Xa 


8 


0.12907 (2) 


0.8648 (1) 


0.81028 (3) 


2.166 (1) 


16 


0.08611 (2) 


2.1984 (3) 


0.76386 (5) 


8.439 (6) 


24 


0.06785 (3) 


3.4141 (6) 


0.73829 (8) 


18.39 (2) 


32 


0.05718 (3) 


4.5578 (8) 


0.72108 (9) 


31.66 (5) 


40 


0.05006 (4) 


5.6466 (11) 


0.70812 (11) 


48.15 (9) 


48 


0.04504 (4) 


6.6993 (14) 


0.69729 (13) 


68.11 (14) 


56 


0.04107 (4) 


7.7199 (23) 


0.68869 (16) 


90.80 (24) 


64 


0.03793 (4) 


8.7142 (25) 


0.68130 (15) 


116.3 (3) 


80 


0.03324 (4) 


10.635 (3) 


0.66899 (19) 


176.6 (6) 


96 


0.02987 (4) 


12.477 (6) 


0.65895 (22) 


249.0 (10) 


112 


0.02723 (5) 


14.273 (6) 


0.65095 (26) 


330.7 (14) 


128 


0.02513 (5) 


16.027 (8) 


0.6441 (3) 


422.4 (23) 



Using the nonlinear Levenberg-Marquardt least-squares algorithm, we fitted for the expo- 
nents and amplitudes according to the finite-size-scaling formulas including correction terms, 
as given in Eqs. (fT3j). (HH), (pUj) and (j221)- Thus we can numerically determine X;, X^, and 
thereby the fractal dimension di of the spanning loop, and the fractal dimension da of interior 
of the spanning loop. Comparing the residual of the fits with the number of degrees of 
freedom, we find satisfactory fits including all system sizes L > 8. We obtain Xi = 0.593 (2) 
from the fit of si{L) and X; = 0.595 (2) from the fit of x{L)- These results are consistent 
with the theoretical value Xi ^ 0.593513601. The fit of Sa{L) yields Xa = 0.080 (1). From 
the fit of XaiL), we obtain 2 - 2Xa = 1.84 (1), or X„ = 0.080 (5). These results agree well 
with the theoretical prediction X^ ^ 0.0801085234. 

In these fits, the correction-to-scaling exponent yi was left free. The fits suggest that 
the exponent of the dominant correction to scaling does not assume the expected value 
2 — X21 = —0.43859, but instead yi = —0.75 ± 0.10. On the other hand, we have also 
performed similar simulations of the 0(1) model at the critical point. The fractions si 
and Sa, and the susceptibility- like quantities XhXa were sampled for several system sizes. 

1 n 



The results are shown in Tables |n] Least-squares fit results agree well with the theoretical 
prediction, as listed in Table IVIl For the 0(1) model, only the data for xi allowed a 
reasonably accurate estimate of the leading correction-to-scaling exponent. In this case, we 
find Hi = —0.628 (7), in a good agreement with the expected value 2 — X21 = —0.625. We 
thus have fixed yi at this value in the other fits. 

Furthermore we performed similar simulations of the 0(n) models with n = a/2, n = -\/3 
and n = 2 at their critical points as given in Ref. 3- The fractions si and Sa, and the 
susceptibility-like quantities XhXa were sampled for several system sizes. The results are 
shown in Tables IIIH IIVI and El Also these quantities appear to depend algebraically on the 
system size L, in agreement with the finite-size scaling equations (HH), (fTSj) . and ((221). 



TABLE II: Numerical data for si{L), xK-^)) Sa{L) and Xa{L) for several system sizes L at the 
critical point of the 0(n = 1) model. 



L 


Sl 


XI 


Sa 


Xa 


8 


0.09126 


(2) 


0.43532 


(8) 


0.88355 (2) 


1.1563 


(6) 


16 


0.06015 


(2) 


1.1859 


(2) 


0.85090 (5) 


4.580 


(5) 


24 


0.04678 


(2) 


1.8531 


(3) 


0.83257 (7) 


10.108 (16) 


32 


0.03911 


(2) 


2.4658 


(4) 


0.81987 (9) 


17.682 


(39) 


40 


0.03405 


(3) 


3.0413 


(6) 


0.81015 (12) 


27.27 


(8) 


48 


0.03040 


(3) 


3.5901 


(8) 


0.80234 (14) 


38.73 


(13) 


56 


0.02754 


(3) 


4.1149 


(12) 


0.79617 (15) 


51.68 


(19) 


64 


0.02539 


(4) 


4.6235 


(14) 


0.79038 (20) 


66.98 


(32) 


80 


0.02214 


(5) 


5.5903 


(24) 


0.78077 (27) 


104.0 


(7) 


96 


0.01968 


(5) 


6.5121 


(33) 


0.77384 (28) 


145.3 


(11) 


112 


0.01781 


(6) 


7.4027 


(42) 


0.76799 (35) 


193.0 


(18) 


128 


0.01645 


(6) 


8.2442 


(43) 


0.76213 (42) 


253.2 


(26) 



Using finite size scaling and the Levenberg-Marquardt algorithm, we determined the 
exponents Xi and Xa- The results are summarized in Table IVTl 
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TABLE III: Numerical data for si{L), xi{L), Sa{L) and Xa{L) for several system sizes L at the 
critical point of the 0(n = ^/2) model. 





Sl 


XI 


Sa 


Xa 
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40 


0.04693 


(4) 


5.1183 


(11) 


0.7264 


(1) 


44.12 


(9) 


48 


0.04206 


(4) 


6.0630 


(16) 


0.7164 


(2) 


62.14 


(18) 


56 


0.03835 


(5) 


6.9786 


(18) 


0.7081 


(2) 


82.86 


(22) 


64 


0.03536 


(4) 


7.8623 


(22) 


0.7011 


(2) 


106.5 


(3) 


80 


0.03096 


(5) 


9.580 


(4) 


0.6890 


(2) 


162.4 


(6) 


96 


0.02769 


(5) 


11.222 


(5) 


0.6800 


(3) 


227.0 


(12) 


112 


0.02532 


(6) 


12.825 


(7) 


0.6716 


(3) 


304.4 


(18) 


128 


0.02338 


(6) 


14.370 


(9) 


0.6649 


(4) 


390.7 


(25) 



V. DISCUSSION 

The time-consuming character of simulations of loop models is due to the nonlocal char- 
acter of the loops. We have reduced this problem by splitting the loop weight n in two 
parts, namely n — 1 and 1. Proper summation over both contributions is done by assigning 
color variables to the loops; a sum on all color variables is included in the partition sum. 
The algorithm treats these color variables as dynamical variables which are updated by the 
Monte Carlo process. The idea to use an additional color variable for each loop was already 
used in a context unrelated to Monte Carlo methods, e.g. in Ref. 13- The presence of loops 
of weight 1 in such a configuration then leads, at least locally, to a system that behaves 
precisely as an Ising configuration. Thus, the system may be locally updated by means of 
Metropolis-like Monte Carlo steps. Care should be taken to leave the loops of weight n — 1 
unchanged, because it would violate the condition of detailed balance. The procedure given 
in Sec. IIIII satisfies this constraint. 

1 o 



TABLE IV: Numerical data for si{L), xii^), Sa{L) and Xa{L) of the critical 0(n = \/3) model for 
several system sizes L. 
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(12) 
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41.16 


(5) 


A n 


0.06111 


(4) 


7.4130 


(17) 


0.6522 


(1) 


62.44 


(9) 


48 


0.05500 


(4) 


8.8463 


(23) 


0.6408 


(1) 


87.2 1 


(13) 


56 


0.05035 


(4) 


10.241 


(3) 


0.6314 


(1) 


115.7 


(2) 


64 


0.04665 


(4) 


11.611 


(4) 


0.6233 


(1) 


148.0 


(3) 


80 


0.04112 


(5) 


14.272 


(6) 


0.6098 


(2) 


223.7 


(6) 


96 


0.03701 


(4) 


16.885 


(7) 


0.5993 


(2) 


311.0 


(9) 


112 


0.03387 


(4) 


19.407 


(9) 


0.5906 


(2) 


411.4 


(12) 


128 


0.03153 


(5) 


21.865 


(11) 


0.5824 


(2) 


531.2 


(19) 



The development of this coloring algorithm was motivated by the possibility to further 
explore the physics of Qin) models. In the course of this work, we realized that it should be 
possible to construct an even more efficient algorithm of a 'cluster' nature. Algorithms of 
th,s type will be presented elsewhere Q. The .co.otmg- tdcU is only useful for „ > 1, For 
< n < 1 the existing algorithm [1^ is, although relatively inefficient, still applicable. 

For the interpretation of the simulation results, it is relevant that we are using periodic 
boundary conditions, and that the mapping between loop and Ising configurations imposes a 
condition of 'evenness' on the loop configurations: a path spanning the periodic boundaries 
must have an even number of intersections with a loop. Therefore, the statistical ensemble 
generated by the algorithm does not coincide with that of Eq. (j3]). The difference is related 
to the boundary conditions and is expected to modify the finite-size behavior, but should 
vanish in the thermodynamic limit. 

The present work is restricted to the 'even' loop configurations. It is, however, possible 
to simulate 'odd' loop configurations by introducing a 'seam' on the dual lattice, a vertical 



TABLE V: Numerical data for s/(-L), xK-^)) Sa{L) and Xa{L) of the critical 0(n = 2) model for 
several system sizes L. 
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(6) 
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(2) 
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(4) 


0.50468 


(5) 


232.70 


(14) 


80 
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(3) 


26.676 
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0.49057 
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345.66 


(30) 


96 
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(2) 
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(6) 
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(37) 
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(3) 


37.828 


(11) 


0.47013 
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(5) 
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0.05410 


(3) 


43.392 


(13) 


0.46223 


(7) 
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(7) 
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0.04418 


(3) 


65.680 


(26) 


0.43916 (10) 
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(2) 



TABLE VL Results for the exponents Xi and Xa for five 0(n) models, as obtained from least- 
squares fits as described in the text. 
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from Xa 
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0.625 
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0.73 


(2) 


0.0518 


(3) 
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(1) 
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(1) 


3/5 


0.796 


(5) 
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(1) 
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(4) 
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(1) 
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1.84 


(1) 


V3 


0.571 


(1) 
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(1) 


0.095 


(1) 


0.0952 


1.81 


(1) 


2 


0.4997 


(3) 


1/2 
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(5) 


0.1249 


(4) 


1/8 


1.747 


(2) 



column of horizontal antiferromagnetic Ising bonds spanning the system. For these anti- 
ferromagnetic bonds we use the rule that there is a loop segment if and only if the two 
associated dual spins are of the same sign. Horizontal and vertical seams can be introduced 
independently, as prescribed by the class of loop configurations that is to be sampled. 
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FIG. 2: Distribution Pi of loops of length / on logarithmic scales, for the critical 0(n) model with 
n = 1.5 and size L = 128. The dashed line shows a power law decay with exponent -2.42198318, 
which is the theoretical asymptotic value for the infinite system. 
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FIG. 3: Fraction si of lattice edges covered by the spanning loop, versus system size L for the 
critical 0{n) model with n = 1.5, on logarithmic scales. The curve is added as a guide to the eye, 
and estimated error bars are smaller than the size of the symbols. 
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FIG. 4: Susceptibility-like quantity xi versus system size L for the critical 0(n) model with n = 1.5 
on logarithmic scales. The curve is added as a guide to the eye, and estimated error bars are smaller 
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